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We report on the results ol the search for extremely-high energy (EHE) neutrinos with energies 
above 10 7 GeV obtained with the partially (~30%) constructed IceCube in 2007. From the absence 
of signal events in the sample of 242.1 days of effective livetime, we derive a 90% C.L. model 
independent differential upper limit based on the number of signal events per energy decade at 
E 2 (f> Vl! + Vfi + l , T ~ 1.4 x 10 -6 GeV cm -2 sec -1 sr _1 for neutrinos in the energy range from 3 x 10 r to 
3 x 10 9 GeV. 

PACS numbers: 98.70.Sa, 95.85.Ry 



I. INTRODUCTION 



Detection of extremely-high energy (EHE) neutrinos 
with energies greater than 10 7 GeV may shed light on 
the longstanding puzzle of the origin of EHE cosmic- 
rays [l|, l2[ . Several observational results have indicated 
that these EHE cosmic rays (EHECRs) are of extra- 
galactic origin 3]. Further elucidation of their produc- 
tion mechanism by EHECR observation is, however, lim- 
ited because the collisions of EHECR with the cosmic 
microwave background photons, known as the Greisen- 
Zatsepin-Kuzmin (GZK) mechanism [4], prevent EHE- 
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CRs from propagating over cosmological distances with- 
out losing a sizable fraction of their energy. On the 
other hand, cosmogenic neutrinos Q produced by the 
GZK mechanism via photo-produced n meson decay as 
7r ± — > iF^v^ — > e^UeU^ carry information on the EHECR 
source evolution and the maximum energy of EHECRs 
at their production sites [f|. 

Detection of these EHE neutrinos is an experimental 
challenge because the very low EHE neutrino fluxes re- 
quire a very large detector. The large size of the IceCube 
neutrino observatory , currently under construction at 
the geographic South Pole, will make it more effective 
than previous experiments in the search for these neu- 
trinos [H, Interactions of u^, v e , and v T and their 
antiparticles are observed through the Cherenkov radi- 
ation emitted by secondary particles. In the following, 
we do not distinguish between v and F; the simulations 
and sensitivity calculations assume an equal mixture of 
particles and antiparticles. 

In this paper we will describe the first results of a 
search for signatures of cosmogenic neutrinos in the 2007 
data acquired by the partially constructed IceCube neu- 
trino observatory. This analysis selects events which pro- 
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duce a large amount of light in the detector. Based on 
simple criteria, such as the total number of observed 
Cherenkov photons and the results of reconstruction al- 
gorithms, it selects candidate neutrino events. Although 
Ufi, u e , and v T interactions look very different in IccCubc, 
the selection criteria arc sensitive to all three flavors. 



II. THE ICECUBE DETECTOR 

IceCube is a cubic-kilometer, high-energy cosmic neu- 
trino telescope which is currently under construction. It 
uses the 2800 m thick glacial ice as a Cherenkov medium. 
Cherenkov photons emitted by rclativistic charged par- 
ticles, notably muons, electrons, and taus produced in 
charged current interactions and their secondaries, are 
detected by an array of photon sensors, known as Digital 
Optical Modules (DOMs) [ijj. The DOMs deep below 
the ice surface are deployed along electrical cable bun- 
dles that carry power and communication between the 
DOMs and surface electronics. The cable assemblies, of- 
ten called strings, are lowered into holes drilled to a depth 
of 2450 m. The DOMs, spaced at intervals of 17 m, oc- 
cupy the bottom 1000 meters of each string. The strings 
are arranged in a hexagonal lattice pattern with a spac- 
ing of approximately 125 m. DOMs are also frozen into 
tanks located at the surface near the top of each hole. 
The tanks constitute an air shower array called IceTop 

The DOMs enclose a down-looking 25 cm photomul- 
tiplier tube (PMT) [12J with data acquisition and cali- 
bration electronics, light emitting diodes for calibration, 
and also data compression, communications, and control 
hardware 

.10] 

in a 35 cm diameter pressure sphere. Al- 
most all of the PMTs are run at a gain of 10 7 ; PMT satu- 
ration effects become important at signal levels of about 
5000 photoelectrons in a single DOM in 50 ns. When 
the DOM detects a photoelectron, it initiates an acquisi- 
tion cycle, recording the PMT output with two waveform 
digitizer systems. The first system samples every 3.3 ns 
for 400 ns, with 14 bits of dynamic range. The second 
system samples every 25 ns for 6.4 /is, with 10 bits of 
dynamic range. The data acquisition system is designed 
such that the first system is sensitive to a bright photon 
source at close distance and the second system captures 
signal induced by photons emitted at large distance. This 
analysis uses the total number of photoelectrons detected 
by the PMTs as a measure of the event energy. For each 
DOM, the charge used is the one from whichever system 
recorded a larger number of photoelectrons. Because of 
the significant DOM-to-DOM differences in saturation 
behavior, the current analysis does not attempt to cor- 
rect for PMT saturation. So, the signals from brightly 
illuminated DOMs are naturally truncated. 



III. DATA AND SIMULATION 

This analysis uses data collected from May 2007 
through April 2008, when IceCube consisted of 22 in-ice 
strings (IC-22; 1320 DOMs) and 52 IceTop tanks. In or- 
der to greatly reduce random noise from radioactivity, in 
IC22 the DOMs only recorded signal waveforms when a 
local coincidence (LC) condition was satisfied, i.e. when 
an adjoining or next-to-nearest neighbor DOM was trig- 
gered within ±1 (is. In 2007, the trigger selected time pe- 
riods when 8 or more DOMs recorded LC signals within 
5 /xs; when this happened, all hits within a 20/is window 
were stored as an event. The average trigger rate was 
about 550 Hz. The high-multiplicity event sample used 
in this analysis imposes an additional condition requiring 
NDOM > 80, where NDOM is the number of hit DOMs 
in an event. The average high- multiplicity event rate was 
approximately 1.5 Hz with a seasonal variation of 17%. A 
total of 3.2 x 10 events were tagged as high-multiplicity 
during the effective livetime of 242.1 days (excluding the 
periods of unstable operation) . 

The high multiplicity cut reduces the data by a fac- 
tor of ~ 3 x 10~ 3 while preserving approximately 70% of 
the GZK neutrinos with projected trajectories that pass 
within 880 m of the center of IceCube. Here, and be- 
low, the GZK signal rates are based on the GZK spectra 
and flux calculated by Ref. @ assuming an all-proton 
composition with a moderately strong source evolution, 
(z + l) m with m = 4 extending to z = 4.0. Neutrino oscil- 
lations modify the neutrino flavor ratio over the cosmo- 
logical distances they travel and the fluxes at the Earth 
were calculated as in Ref. [l3[ . Note that the flavor ratio 
v e : : v T of cosmogenic neutrinos at the Earth is dif- 
ferent from 1:1:1 as primary energy spectra of v e and 
produced by the GZK mechanism are different because of 
a significant contribution of v e from neutron decay. This 
enhancement was included in the GZK neutrino flux cal- 
culations used here. 

EHE neutrinos were simulated with the JULIcT pack- 
age [9( to generate and propagate the neutrinos through 
the Earth. All three flavors of neutrinos were simulated 
with energies between 10 5 and 10 11 GeV. The result- 
ing secondary muons and taus produced in the neutrino 
interactions are propagated through the rock and ice 
near the IceCube volume, also by JULIeT. Hadronic and 
electromagnetic showers are also simulated; all of these 
showers are treated as point sources, without account- 
ing for the LPM effect. The background muon bundles 
from cosmic-rays in the energy range 10 6 GeV to 10 10 
GeV were generated using CORSIKA [l4| version 6.720 
with the SIBYLL 2.1 hadronic interaction model or with 
QGSJET II, without charm production |l5]. The un- 
certain prompt muon component from charm decay may 
contribute to the background events [l6[. The muons 
were propagated through the Earth using MMC [TtJ • De- 
partures of observed data distributions from those of the 
CORSIKA based background events prompted us to also 
develop a phenomenological background model based on 
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fits to data. Emission of Cherenkov photons and their 
propagation in the ice was simulated by the Photonics 
package [l8j]. Measurement of the absolute number of 
Cherenkov photons is important in the EHE neutrino 
search as it closely relates to the energy of the high energy 
muons, taus or electrons produced by EHE neutrinos. 
Therefore the detection efficiency of the DOMs must be 
understood with good precision. The primary element, 
the PMT, is calibrated in the laboratory using a nitro- 
gen laser to measure the photon detection efficiency [l2j ■ 
This bare PMT data are used in a simulation package 
which propagates photons inside the glass sphere and the 
optical gel to the photo-cathode surface. The DOM simu- 
lation is followed by waveform calibration and the trigger 
condition is included in the simulation chain. 




6 6.5 7 7.5 8 8.5 9 9.5 10 10.5 11 
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IV. ANALYSIS 

A. Extremely high energy event signatures and the 
initial event filter 

The event signatures from v^, v e , and v T are very dif- 
ferent. IceCube mainly detects cosmogenic neutrinos by 
the signals from the secondary muons and taus generated 
in the neutrino interaction in the rock or ice. At high en- 
ergies, these particles are seen in the detector as series 
of energetic cascades from radiative energy loss processes 
such as pair creation, bremsstrahlung, and photonuclear 
interactions, rather than minimum-ionizing tracks. The 
radiative energy losses are approximately proportional 
to the energy of the muon or the tau, and so is the 
Cherenkov light yield. Electron neutrinos produce elec- 
tromagnetic and hadronic showers, which are relatively 
compact sources of Cherenkov light. Muon and tau neu- 
trinos within IceCube will also produce a hadronic shower 
from the struck nucleon, in addition to the muon or tau 
secondary. 

Shown in Fig. [T] is the simulated distribution of the to- 
tal number of photo-electrons per event (NPE) recorded 
by the the IC-22 detector as a function of the simulated 
true muon energy. A clear correlation between NPE 
and the energy of particles measured near IceCube is 
observed. The energies are sampled at a radius of 880 
m from the IceCube center. This definition of energy 
is labeled "in-ice energy" and used throughout this pa- 
per. The visible departure from linearity for large NPE 
stems from the saturation of the detector during signal 
capture. Approximately 30% of EHE signal events are 
due to neutrino interactions inside the IceCube detector 
volume initiating a hadronic or electromagnetic cascade. 
The correlation between NPE and incoming neutrino en- 
ergy also holds for these events. Electron neutrinos are 
detectable via this channel. 

Because the energy spectrum of background atmo- 
spheric muons (both single muons and bundles) falls 
steeply with energy, the GZK neutrino flux should domi- 
nate over background in the high NPE region. Since the 



FIG. 1: Event distribution from Monte Carlo simulations 
of single muons with the IC-22 detector configurations in a 
plane of NPE and simulated true energy. The muon energy is 
given when the muon is 880 meters from the IceCube center 
(in-ice energy). The 80 DOM multiplicity cut (level 1 cut) is 
applied. The charged lepton energy distribution is assumed 
to follow E~ x in this plot for illustrative purposes. Only par- 
ticles with trajectories intersecting within 880 m from the 
center of IceCube array are considered in the plots. More 
distant events do not contribute to the data sample. 



through-going muons and taus induced by EHE neutrinos 
enter into the IceCube volume mainly horizontally [H, |9| , 
the signal search criteria are chosen to favor roughly hor- 
izontal high NPE events. 

The high-multiplicity NDOM > 80 sample is domi- 
nated by atmospheric background muons. The next step 
of the analysis selects events with NPE > 10 4 . This re- 
duces the background by three orders of magnitude, leav- 
ing 6528 events, still dominated by background, while the 
GZK signal reduction is ~24%. 

Table U summarizes the number of events remaining at 
each level of the initial filtering. In order to estimate the 
background in the very high energy region, the simulated 
data are compared to the experimental data in the region 
10 4 < NPE < 10 5 . The present analysis follows the blind 
analysis technique. In keeping with the IceCube blind- 
ness policy, events with NPE > 10 5 were not used for 
determining the background or setting cuts. This NPE 
threshold was chosen so that the possible contribution 
from signal events in the studied sample was negligible. 



B. High-energy muon background 

Bundles of muons generated in cosmic-ray air showers 
are the major background for the EHE neutrino signal 
search, because multiple muon tracks with a small geo- 
metrical separation resemble a single high energy muon 
in the IceCube detector. The multiplicity, energy distri- 
bution, and separation distances for these muon bundles 
are not fully understood. Two independent Monte Carlo 
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TABLE I: Number of events at different filter levels for 242.1 days in 2007. The simulation predictions for the atmospheric 
muon background using the CORS1KA-SIBYLL package, the empirical model, and that for the GZK cosmogenic neutrino 
model are also listed for comparison. Errors shown here are statistical only. Refer to Sec. II V Al for the definitions. 



Filter levels 


observational data empirical model 


CORSIKA (iron) 


CORSIKA (proton) 


signal (GZK1 [6]) 


level 1 (NDOM > 80 
level 2 (NPE > 10 4 ) 


) 3.195 x 10 7 

6528 (6.82 ± 0.42) xlO 3 


(1.84± 0.08) xlO 7 
(1.09 ± 0.09) xlO 4 


(7.71 ± 0.45) x 10 6 
(1.63 ± 0.17) xlO 3 


(886 ± 8.9) xlO -3 
(670 ± 7.5) xl0~ 3 




FIG. 2: Event distributions for NPE, cosine of reconstructed zenith angle, and the NPE-weighted mean depth of event (zcog) 
for observational and the background MC simulation data. The black dots represent observational data after the NPE > 10 4 
cut, red for CORSIKA proton (SIBYLL), magenta for CORSIKA iron (SIBYLL). Green shaded regions represent distributions 
obtained with the empirical model with the size of shade expressing the uncertainty of the model. See text for the details. 



simulations are carried out to estimate the muon-bundle 
background in this EHE neutrino signal search. 



The first is the full cosmic-ray air shower simulation 
with light and heav y io n primaries using the CORSIKA 
(SIBYLL) package d4|. Two extreme cases of compo- 
sition are used to address the event rate variation due 
to the uncertainty in the primary cosmic-ray mass pop- 
ulation. While the full air shower simulation includes 
a calculation of the production spectra of the multiple 
muons from meson decay, the simulation still introduces 
a large uncertainty because both the primary composi- 
tion at relevant cosmic-ray energies (> 10 7 GeV) and the 
hadronic interaction model are highly uncertain. 



The second simulation uses a model relying on a phe- 
nomenological fit to part of the experimental high energy 
data. This empirical model approximates multiple muon 
tracks in an event by a single high energy muon which is 
adequate at very high energies for the variables used in 
this analysis. The single muon approximation predicts 
larger fluctuations in NPE due to radiative energy losses 
of energetic muons, giving a rather conservative estimate 
of the background passing rate. The two independent 
sets of simulations with top-down and bottom-up ap- 
proaches to describe the observational data complement 
each other improving the reliability of the background 
estimation. 




log (Total NPE) z COG [m] 



FIG. 3: Event distributions from the level-3 samples as func- 
tions of NPE (left) and zcoa (right). The black dots represent 
observational data, green boxes represent the empirical model 
including uncertainty. Red and magenta lines are CORSIKA 
samples with SIBYLL interaction model and proton and iron 
primaries, respectively. The left panel also includes the ex- 
pected NPE distribution of events induced by the cosmogenic 
neutrinos @] shown by the blue line for reference. 



1. Background estimation with CORSIKA 

Fig. [5] shows distributions for data and simulations 
at level-2 for NPE, the reconstructed zenith angle (9), 
and the center-of-gravity depth of the events (zcog)- 
The CORSIKA NPE distribution is extrapolated to the 
higher NPE region. Extrapolation was necessary mainly 
because of a lack of simulated CORSIKA events at pri- 
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mary cosmic ray energies above 10 10 GeV. The extrap- 
olation accounts for the observed GZK cutoff at an en- 
ergy around 5 x 10 19 eV. The NPE-weighted LineFit al- 
gorithm was used to reconstruct zenith angle in this ini- 
tial study. The NPE-weighted LineFit is a simple min- 
imization of x 2 — ^iNPEi(fi — rcoG ~ tiil) 2 , where 
ti and NPEi represent, respectively, the time of the 
first photoelectron and the number of photoelectrons 
recorded by the i'th DOM at the position fl and fcoc = 
( XiNPEj xj s.jVPg, y> TjjNPEj Zi \ • NPE- weighted 

\ YliNPEi ' T,iNPEi ' T,iNPEi > h 1 K 1 ^ rEj weignicu 

position of the center-of-gravity of the hits. The fit ig- 
nores the geometry of the Cherenkov cone and the opti- 
cal properties of the medium and assumes light traveling 
with a velocity v along a 1-dimensional path through the 
detector, passing through the center-of-gravity. 

The measured event rates are close to the simulated 
rates based on CORSIK A/SIB YLL with iron primaries 
and above those based on CORSIK A/SIB YLL proton 
data in most regions. A significant discrepancy can 
be found in the rate of events with cos 9 < 0.3, i.e. 
events reconstructed as horizontal or up-going, which is 
largely underestimated. Replacing SIBYLL with other 
hadronization models (e.g. QGSJET-II) does not change 
this behavior. The discrepancy may be due to a combina- 
tion of uncertainties in the hadronic interaction models, 
cosmic-ray flux, and Cherenkov photon propagation in 
the glacial ice. Since the horizon is the key region for the 
EHE neutrino search, the background estimations were 
supplemented by an empirical model fit to a subsample 
of the data. 



2. Construction of an empirical model 

The empirical model is optimized to match the level- 
2 experimental data (10 4 < NPE < 10 5 ). The possible 
signal region (NPE > 10 5 ) is not used to avoid bias. The 
model provides a relation between the NPE of an event 
and the cosmic-ray primary energy. Its convolution with 
the cosmic ray flux then gives the event rate with a given 
NPE. The cosmic-ray flux used in the present analysis is 
taken from the compilation in Ref. 0. 

The model is based on the so-called Elbert formula [l9| 
which parametrizes the mean multiplicity of muons with 
energies above a certain threshold E^: 

E cosfl' \ E J \ E Q J K ' 
E T = 14.5 GeV 

where A, Eq, and 9' are the mass number, the energy, 
and the zenith angle of the primary cosmic-ray [20j . The 
energy weighted integration of the formula relates the 
total energy carried by a muon bundle to the primary 
cosmic-ray energy Eq, 

E /A 

E^ ee I — —Ei,dEu 



assuming AE^/Eq <C 1. Here, e is empirically deter- 
mined by fit to the observed data. Assuming its cor- 
responding energy at the IceCube depth, e ln - lce 5 i s in- 
dependent of zenith angle, e (and thereby E^) can be 
calculated as a function of zenith angle by taking into 
account the energy loss during propagation through the 
Earth. The optimization of the two parameters a and 
£ m-ice j g p er f ormec } by comparing the observed data to 
simulation of a single high energy muon with energy of 
in the NPE and zenith angle space, independently. 
A = 1 is assumed in the optimization. The event distri- 
butions derived from the empirical model with optimized 
parameters (a = 1.97 and e ln - lce = 1500 GeV) are given 
in Fig. [2J The green shaded region in the plot is ob- 
tained by allowing the model parameters to vary within 
± 1 a from their optimized values. The discrepancies 
of zcog at large depths for the empirical model and at 
small depths for CORSIKA/iron seem to be due to verti- 
cal, down-going events because a restriction of the zenith 
angle, cos 9 < 0.8, improves the agreement in both cases. 
Since the majority of the EHE neutrino induced events is 
close to the horizon we can discard all events with cos 9 
< 0.8 without significant loss of signal efficiency (level-3 
cut). The resulting distributions are shown in Fig. El 

C. Search for EHE cosmogenic neutrino signal 

The level-4 cut to eliminate the muon background is 
carried out in the NPE-cos 9 (NZ) plane. In accordance 
with the requirements of blindness, the cuts are finalized 
on simulated events alone without referring to real data. 
Because the optical properties of the glacial ice vary sig- 
nificantly with depth [2l| . and because the changing ab- 
sorption and scattering lengths affect what IceCube ob- 
serves, the final cuts are chosen to be depth dependent. 
The cuts are chosen based on the depth of the weighted 
center of gravity of the event, zqog- The distribution 
of events in the NZ plane depends on zcog ■ We divide 
the events into two groups according to their zcog a s 
follows: 

region A: —250 < zqog < —50 m and zcog > 50 m 
region B: zqog < —250 m and —50 < zqog < 50 m 

As seen in Fig. region B contains a large number 
of horizontal and up-going mis- reconstructed background 
events, whereas the fraction of such events in region A 
is very small. Fig. [4] also shows the distributions of the 
experimental data and the simulated GZK neutrino in- 
duced signal events. The latter clearly accumulate near 
the horizontal direction regardless of the zqog position 
and have on average larger NPE than the background 
sample. The selection criteria to separate signal from 
background are determined for region A and B sepa- 
rately: For each bin of cos 9 (width 0.1), a threshold 
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TABLE II: Number of events at analysis filter levels for 242.1 days in 2007. The simulation predictions for the atmospheric 
muon background using the CORSIKA-SIBYLL package, the empirical model, and that for the GZK cosmogenic neutrino 
model are also listed for comparison. Errors shown here are statistical only. See Sec. II V Bl and IIV CI for details. 



Analysis filter levels 


observational data empirical model 


CORSIKA (iron) 


CORSIKA (proton) 


signal (GZKl [6]) 


level 3 (cos(0) < 0.8) 
level 4 (EHE v search) 


2014 (2.65 ± 0.21) xlO 3 
(6.32 ± 1.37) xlO" 4 


(2.68 ± 0.19) xlO 3 
(4.18 ± 1.29) xlO -4 


(4.16 ± 0.40) xlO 2 
(1.44± 0.58) xlO -4 


(620 ± 7.3) xlO -3 
(155 ± 1.4)xl0 -3 




log M (NPE) log io (NPE) logJNPE) log^NPE) 

FIG. 4: Event number distributions passing the level-2 selection cut (NPE> 10 4 ) of the experimental data (left), the background 
from the empirical model (middle left), the background from CORSIKA SIBYLL with iron primaries (middle right), and the 
signal (right) on the NZ plane at the IceCube depth. The upper (lower) panels show the distributions in the region A (B). The 
GZK neutrino flux [6J determines the event intensity in the signal MC plot, adding all three flavors of neutrinos. The series of 
thick lines in each panel indicates the level-3 (cos(#) < 0.8) and the final level-4 cuts. 



TABLE III: Expected event numbers passing the final level-4 
selection criteria in the 2007 IC-22 observation. Models in- 
clude the GZK models [HHm and the Z-burst model 
The predictions are normalized to a livetime of 242.1 days. 
Signal event numbers represent the sum over all three neutrino 
flavors. The first uncertainty is the statistical uncertainty de- 
termined by signal simulation statistics, and the second is the 
total systematic uncertainty from sources discussed in Sec.lVl 



Models 


Number of Events per 242.1 days 


GZKl [6J 
GZK2 [22] 
GZK3 [23] 
Z-burst [29] 


(155 ± 1.4 tit )xl0 -3 
(248 ± 2.3 til )xl0" 3 
(83 ± 0.8 +2 3 )xl0~ 3 
(398 ± 3.4 ttl )xl0" 3 



NPE is set such that the number of background events 
above the threshold is less than 10~ 4 . Tighter cuts to 
further reduce the background would also reduce the sig- 
nal to an undesirable degree. The NPE thresholds in all 
the zenith angle bins are then connected to each other to 
form a series of lines on the NZ plane, defining the final 
level-4 cut, as drawn in Fig. 3J The cuts were optimized 



using the empirical model for background simulations, so 
we used extrapolated CORSIK A/iron data as a check of 
the final background level. Figure 2] also shows the dis- 
tribution of background events in the NZ plane from the 
extrapolation. Table [TT1 summarizes the number of events 
remaining in the analysis after each of the cut levels. 

The effective area as a function of energy at the Earth 
surface for each neutrino flavor is shown in the left panel 
of Fig. [5j averaged over all solid angles. The area in- 
creases with the energy owing to the increasing neutrino 
interaction cross sections and the increased probability 
of observing the interactions, which is different for each 
flavor. 

At low energies, most of the v r signal comes from 
events where the v T interacts in the detector, or the r 
decays in it. At higher energies, r energy loss becomes 
large enough that through-going taus also pass the cuts. 
Contributions from and v T dominate over v e in the 
energy range above ~ 10 8 GeV, as the secondary pro- 
duced muons and taus can travel long distances to reach 
the detection volume. This trend is reversed at lower en- 
ergy where tau and muon energy losses are smaller, and 
j/ e 's can deposit all of their energy into the detector vol- 
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6.5 7 7.5 8 8.5 9 9.5 1010.5 
log, (Energy/GeV) 



6 6.5 7 7.5 8 8.5 9 9.5 1010.5 
log (in-ice energy/GeV) 



FIG. 5: The effective area of IC-22 for EHE neutrino search. 
The left panel shows the 4-7T solid angle averaged area as a 
function of neutrino energy at the Earth surface. The right 
panel shows the corresponding effective area for particles at 
880 m from the IceCube center entering into the IC-22 fidu- 
cial volume. Muons and taus in this plot are secondary par- 
ticles produced by neutrinos before reaching the neighbor- 
hood of the detector array. The energy here are defined as 
in-ice energy. 



TABLE IV: List of the statistical and systematic errors. The 
signal rate is estimated by assuming the high evolution flux 
(m,Z max ) = (4,4) in Ref @. 



Error source 


background 
rate 


signal (GZK) 
rate 


Sltaf i^ripfil prrnr 


±22% 


±0.9% 


Detector sensitivity 




±8% 


Yearly variation 


±17% 




Empirical model 


±99/-59% 




Background model dep. 


±15 % 




NPE yield 




±0/-21% 


Neutrino cross section 




±9% 


Photo-nuclear interaction 




±10% 


LPM effect 




±1% 


Total 


±22% (stat.) 
±102/-63% (sys.) 


±0.9% (stat.) 
±16/-26% (sys.) 



ume. The effective area for v T is larger than that for 
at low energies, because of the events where taus decay 
inside the detector. At the highest energies, because of 
the larger mass of taus and increase of the tau decay time 
with energy, the tau range is longer than that of muons, 
leading to a larger effective area. 

The right panel in Fig. [5] shows the effective area as a 
function of the in-ice energy ("in-ice area"). It represents 
the probability of detection of incoming particles with 
the present analysis. The area for incident muons and 
taus gradually increases with energy but is limited essen- 
tially by the physical cross section of the IC-22 array, ~ 
0.5 km 2 . Because the Cherenkov yield of taus is smaller 
than muons with the same energy due to the smaller ra- 
diative energy loss, the detection probability of incident 
taus is lower, leading to the smaller in-ice area. Incom- 
ing neutrinos must interact to yield Cherenkov light to be 
detected. Therefore the neutrino effective area becomes 
much smaller than that for muons or taus. 

The expected number of signal events for various neu- 
trino production models after the level-4 cut are summa- 
rized in TABLE Mil GZK1 [6] represents the case of a 
moderately strong source evolution, (z+ l) m with m = 4 
extending to z = 4.0, while GZK2 [22| assumes m = 5 
up to z — 2.0, and GZK3 [23| uses m = 3 with a slightly 
different parametrization and a cut-off structure. 



V. THE SYSTEMATICS 

This search is based on the event-wise NPE and recon- 
structed zenith angle. The main systematic uncertainties 
derive from 1) the necessity to extrapolate the empirical 
fit to data by approximately an order of magnitude in 
NPE to estimate the background rate at the highest en- 
ergies and from 2) the uncertainty of the absolute NPE 



scale. Table IIVI lists the sources of statistical and sys- 
tematics errors. 



A. Uncertainties in the background rate estimation 

The largest uncertainty in the background rate esti- 
mate arises from the fact that the parameters of the 
empirical model were optimized for the observed events 
with 10 4 < NPE < 10 5 after level 2 selection. The lim- 
ited statistics of this sample results in uncertainties on 
the parameters. The model was then extrapolated to a 
higher NPE region for the determination of the level-4 
cut. 

Allowing the parameters to vary within ± 1 a changes 
the background rate by between —59% (for the softest 
possible NPE spectrum after the level-4 cuts) and ± 99% 
(for the hardest possible NPE spectrum). Uncertain- 
ties in the detector sensitivity are incorporated by the 
parametrization. The difference in the background level 
estimated with the extrapolated CORSIKA/iron and the 
empirical model can be taken to indicate the level of sys- 
tematic uncertainty due to model dependence. This un- 
certainty is approximately ±15% and can be assumed to 
include the possible contribution from charm decay. An 
uncertainty associated with the high energy hadronic in- 
teraction model is evaluated using simulated muon bun- 
dle intensity from SIBYLL and QGSJET II with iron 
primaries and found to be ± 4%, which is negligible. An 
additional uncertainty of ~17% arises from the seasonal 
variation of the atmospheric muon rate as the signal se- 
lection criteria are based upon the season-averaged data. 
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B. Uncertainties of the signal rate estimate 

The uncertainty in the relationship between measured 
NPE and the energies of charged particles is the largest 
systematic error affecting the signal event rate. It is the 
consequence of our limited understanding of the detec- 
tor sensitivity, the photon propagation in ice, and the 
detector response to bright signals. It is evaluated using 
absolutely calibrated in situ light sources and amounts to 
a possible overestimation of NPE in simulation by 18.5%, 
which leads to decrease of the GZK signal rate by ~21%. 

Uncertainties of the relevant particle interactions in 
the EHE regime also add systematics in the signal rate 
estimation. The expected event rate scales nearly lin- 
early with the neutrino-nucleon inelastic cross section in 
the EHE range. This scaling has been confirmed by 
numerical studies, artificially increasing the cross sec- 
tion. The cross section uncertainty has been recently 
reduced to be around ±9% with the inclusion of the 
most recent data from HERA and modern parton dis- 
tribution functions 24]. Another systematic error arises 
from the photonuclear cross section of EHE muons and 
taus. The present calculation used the model by Bugaev 
and Shlepin (25|, rewritten in Ref. (26|, that includes a 
relatively reliable soft nonperturbative component and a 
less certain hard perturbative part. Ignoring the hard 
component in the simulation gives the most conservative 
estimate of the uncertainty and leads to a 10% event 
rate increase. The suppression of bremsstrahlung and 
pair production due to the LPM effect (27|, for the rele- 
vant electron energies of 10 9-10 GeV, increases the effec- 
tive radiation length of the electromagnetic cascade to 
O(30 ~ 100) m from ~36 cm (28|. Because the value 
is still comparable to the IceCube DOM separation, and 
the contribution from v e constitutes <20% of the total 
event rate in this energy range, the LPM effect has a neg- 
ligible impact on the event rate. This has been confirmed 
by a special simulation study on v e including the LPM 
cascade elongation. 



VI. RESULTS 

No events are observed in the final data sample taken 
in 2007 with a livetime of 242.1 days when applying the 
final level-4 selection criteria, which is consistent with the 
expected number of background events of 6.3 x 10 -4 . We 
choose to present the resulting all flavor EHE neutrino 
upper limit in the quasi-differential form independent of 
the neutrino production model. Assuming full mixing 
due to oscillations, the experimental 90% confidence level 
upper limit is obtained by setting 2.44 events 34] for an 
upper bound of the number of events observed with bin 
width of a decade of energy with condition that energy 
dependence of neutrino flux multiplied by the effective 
area behaves as ~ 1/E 31]. This limit is presented in 
Fig. |U including the systematic errors. The plot indicates 
that the EHE neutrino search by the IceCube observatory 
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FIG. 6: The all flavor neutrino flux differential limit from 
the IC-22 EHE analysis (filled circles) . The systematic errors 
are included. Also the various model predictions are shown 
for comparison: GZK model 1 .6] (short dashed line), GZK 
model 2 [H (dotted line), GZK model 3 HI (long dashed 
line), Z-burst model [2!| (dashed dot line). The model inde- 
pendent differential upper limits by other experiments are also 
shown for Auger [§3] (open triangles), RICE [3lJ (crosses), 
ANITA [H (open squares), AMANDA [jjjj (rhombi). Limits 
from other experiments are converted to the all flavor limit 
assuming full mixing neutrino oscillations and 90% C.L when 
necessary. 



is most sensitive to the neutrinos with energies on the 
Earth surface ranging between about 10 s and 10 9 GeV. 
The absence of signal events in the sample of 242.1 days of 
effective livetime results in a 90% C.L. differential upper 
limit on the neutrino flux of E 2 <f> Uc+u ^ +UT ~ 1.4 x 10 -6 
GeV cm -2 sec^ 1 sr _1 on average for neutrinos with an 
energy of 3 x 10 7 < E < 3 x 10 9 GeV. Here (t>„ e+Vli +u T 
denotes the differential flux of the sum over all three neu- 
trino flavors, i.e. number of neutrinos per unit energy, 
area, time and solid angle. 

The quasi-differential limit in Fig. [6] takes into account 
the systematic uncertainties. The background rate stays 
negligible O(10~ 3 ) even including the systematic uncer- 
tainty and the resultant upper limit is unchanged. The 
signal rate uncertainty is strongly dominated by the un- 
certainty of the NPE yield which influences the number 
of expected signal events as a function of the neutrino 
energy. The upper limit is calculated by reducing NPE 
by 18.5% in the signal simulation to account for this fac- 
tor. All the other sources of systematic error only slightly 
change the signal passing rate; they are independent of 
energy. They are included in the analysis by uniformly 
scaling the effective area in the limit calculation. 

The present limit is approximately a factor of 20-30 
higher than the intensity range expected in the GZK cos- 
mogenic neutrino production models 0, [H, , as one 
can see in Fig. [51 The current limit for 242.1 days of ob- 
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servation is comparable to the Auger [3(| and HiRes [35[ 
bounds by their multiple year operation. 

VII. CONCLUSIONS 

The present work has demonstrated that the Ice- 
Cube neutrino observatory is capable of searching for 
signatures of EHE cosmogenic neutrinos with relatively 
straightforward event selection methods. The model in- 
dependent differential upper limit obtained with 242.1 
days of observation in 2007, with approximately one 
quarter of the completed detector is E 2 <j) Ue+Ufi+UT ~ 
1.4 x 10 -6 GeV cm~ 2 sec -1 sr^" 1 for neutrinos with an 
energy of 3 x 10 7 < E < 3 x 10 9 GeV. This is approx- 
imately a factor of 20 higher than the predicted GZK 
neutrino flux from relatively strongly evolved sources. In 
the future, data taken by IceCube with 40 to 86 strings 
operating should lead to a detection of cosmogenic neu- 
trinos or a greatly improved limit. 
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